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Abstract 

The study deals with the Payne effect (a substantial decrease in the storage modu- 
lus of a particle-reinforced elastomer with an increase in the amplitude of mechanical 
oscillations). The influence of temperature, concentration of filler and amplitude and 
frequency of strains is analyzed on the mechanical response of filled rubbery poly- 
mers. Constitutive equations are derived using the concept of two interpenetrating 
networks: one comprises semiflexible polymeric chains connected to temporary junc- 
tions, whereas the other is formed by aggregated filler clusters. Adjustable parame- 
ters are found by fitting experimental data for natural rubber, bromobutyl rubber and 
styrene-butadiene rubber reinforced by carbon black and polymeric particles. The 
critical concentration of particles is determined that characterizes transition from an 
ensemble of disjoint clusters to the network of filler. The volume fraction of filler 
corresponding to this transition is found to be close to theoretical predictions based 
on the percolation theory, as well as to experimental data for isolator-conductor 
transition. 



1. Introduction 

This paper is concerned with the Payne effectEhi: a noticeable decrease in the storage mod- 
ulus of filled rubbery polymers with an increase in the amplitude of small-strain oscillations 
in dynamic mechanical tests. Changes in dynamic moduli of filled elastomers are tradi- 
tionally associated with aggregation of filler particles into clusters and networkJl^i. This 
phenomenon has attracted substantial attention of experimentalists in the past decadeH^ 
because of its importance for industrial applications (airplane and automotive tires). 

The strain-dependence of dynamic moduli is conventionally described in terms of the 
Kraus equations^ 

r,, r,/ A£' „ „ A£"(A/A Q ) W 

~ + 1 + (A/A ) 2m ' ~ 1 + (A/A ) 2m ' 1 j 

where E' and E" are storage and loss moduli in a dynamic test with the strain- amplitude A, 
E' Q and Eq are dynamic moduli in a test with an infinitesimally small amplitude of strains, 
and A and m are adjustable parameters. Equations (1) were developed with reference 
to the deagglomeration-reagglomeration mechanism for aggregation of filler clusters. The 
exponent m in Eqs. (1) is fairly well predicted in terms of the connectivity parameter for 
fractal aggregates!. 

Formulas (1) provide an acceptable agreement with experimental data in dynamic tests 
with a fixed frequency of oscillationJIll. However, three shortcomings of the Kraus equations 
may be mentioned: 

1. fitting observations for storage and loss moduli separately results in different values of 
adjustable parametersllS; 

2. Eqs. (1) do not describe changes in dynamic moduli of rubbery polymers with tem- 
perature, volume fraction of filler and frequency of oscillations; 

3. the deagglomeration-reagglomeration concept does not distinguish between the re- 
sponse of an ensemble of disjoint clusters corresponding to small degrees of loading 



and that for a filler network (an agglomerate superstructure^) observed in elastomers 
with high concentrations of reinforcing particles. 

The neglect of the clusters-network transition implies that parameters of Eq. (1) smoothly 
change with the filler concentration. On the other hand, observations in dc conductivity 
tests reveal that conductivity of reinforced rubbery polymers increases by several orders of 
magnitude when the volume fraction of filler reaches some percolation thresholdlih§. This 
contradiction implies the following questions: (i) whether an analog of the isolator-conductor 
transition is observed in dynamic mechanical tests, and (if so), (ii) whether the threshold 
concentrations of filler found in mechanical and dc conductivity experiments are close to one 
another. 

We aim to develop a model which can correctly describe the effects of amplitude and 
frequency of oscillations on dynamic moduli. For this purpose we apply the theory of 
temporary polymeric networksfHEI to characterize the time-dependent behavior of the host 
material and the concept of boundBB and occluded^ rubber to predict interactions between 
the elastomer and filler particles. 

According to conventional theories of transient networks, a rubbery polymer is thought 
of as an incompressible network of long chains connected to temporary junctions. A chain 
whose ends are connected to separate junctions is treated as an active one. Snapping of an 
end of a chain from a junction is modeled as its breakage (transition from its active state to 
the dangling state). When a dangling chain captures a junction, a new active chain arises. 
Breakage and reformation of active chains (occurring at random times) are thought of as 
thermally activated processes. 

Unlike previous studies where all chains are treated as composed of an identical number 
of strands, we suppose that chains contain various numbers of strands, n, and determine 
the distribution of chains by the probability density, p(n), which is assumed to change with 
temperature, the volume fraction of filler and the strain intensity. To reduce the number of 
adjustable parameters in the governing equations, we adopt a discrete analog of the Gaussian 
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distribution, 

p{n) = Aexp[- K ^ 2 ; j (n = l,...,iVi), p(n) = (n = JVj + 1, . . .), (2) 

where A" is the average number of strands in a chain, £ is an analog of the standard deviation 
for the number of strands, N\ is the maximal number of strands in a chain, and the constant 
A is determined by the equality 

Ni 

J^Pin) = 1. 

n=l 

The advantage of Eq. (2) is that given a priori quantities iV and Ai, the distribution of 
chains is determined by the only parameter £. 

Another refinement of the conventional theory of transient networks is that the rates of 
breakage and reformation for active chains are assumed to depend on the number of strands, 
n, which allows relaxation spectra of filled elastomers to be adequately predicted. 

The theories of bound and occluded rubber presume that macro-strains in the host ma- 
terial are constrained by surrounding aggregates of filler. The theory of bound rubber, see 
Refs.@@@ and the references therein, ascribes these constrains to adsorption of macro- 
molecules on the filler surfaces, which results in formation of shells of the host material 
(with the thickness of about 10 nm@) around filler particles and their clusters. Physical 
properties of the surface layers are assumed to substantially differ from those for the unfilled 
elastomer^. In particular, elastic moduli and the resistance to the desorptive action of a sol- 
vent for the shells exceed those for the bulk polymer by several orders of magnitudeEl. NRM 
spectroscopy reveals an internal structure of the shells which consist of an inner domain of 
tightly bound macromolecules and an outer region of loosely bound chains^. 

According to the concept of occluded rubberi, aggregation of isolated clusters of particles 
into a network results in separation of some meso-regions of the host material (occluded 
domains) from the bulk elastomer. When the filler network is relatively rigid, this implies 
shielding of constrained rubbery regions by the network which does not transmit stresses 
to the occluded regions. This results in strong inhomogeneity of macro-deformation. To 



simplify the analysis, the host material is assumed to be split into meso-domains, in a part 
of which micro-strains coincide with the macro-strains in the sample, whereas in the other 
part micro-strains vanish. 

The following scenario is proposed for deformation of a particle-reinforced elastomer. The 
distribution of chains, p(n), in an unfilled rubbery polymer is characterized by a relatively 
small standard deviation E that determines its relaxation spectrum. For very small degrees 
of loading, 0, the parameter E increases with the growth of the filler concentration which 
reflects creation of thin shells containing regions of loosely bound chains (the number of 
strands in chains belonging to these domains is assumed to differ from the "average" number 
for an unbounded rubber). For low concentrations of filler, when clusters do not aggregate 
into a network [the regime of disjoint clusters (DC)], changes in the viscoelastic response of a 
reinforced elastomer are associated with transformations occurring in surface layers around 
the filler particles and their flocks. We suppose that 

1. given a temperature, T, and an amplitude of oscillations, A, the quantity E decreases 
with the degree of loading because of aggregation of isolated particles into clusters, 
whose surface area is substantially less than the sum of surface areas of separated 
particles; 

2. given a degree of loading, 0, and a strain amplitude, A, the standard deviation of the 
number of strands, E, increases with temperature T because of a partial destruction 
of glassy-like shells of tightly linked chains around the filler particles and their clusters 
and transformation of tightly bound layers into loosely bound domains; 

3. given a temperature, T, and a degree of loading, 0, the parameter E weakly decreases 
with the amplitude of oscillations A because of mechanically induced breakage of long 
chains in the regions of loosely bound rubber and their subsequent amalgamation with 
the host elastomer. 

At high concentrations of filler, when clusters of reinforcing particles are aggregated into a 
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network [the regime of network rupture (NR)], changes in the time-dependent response of a 
rubbery polymer are mainly determined by interactions between the host material and the 
network. We suppose that 

1. given a temperature, T, and an amplitude of oscillations, A, the quantity £ strongly 
decreases with the degree of loading because of coalescence of flocks into a fractal 
network, which causes a noticeable reduction in the lateral surface of clusters to which 
regions of loosely bound rubber are adjusted; 

2. given a degree of loading, 0, and a strain amplitude, A, the standard deviation of the 
number of strands, S, decreases with temperature T because of the growth of rate of 
breakage for long chains in the regions of loosely bound rubber and their subsequent 
merging with the bulk polymer; 

3. given a temperature, T, and a degree of loading, 0, the quantity E dramatically 
decreases with the amplitude of oscillations A because of mechanically induced rupture 
of the filler network and its crumbling into separated clusters. 

The objective of this study is to validate this scenario by comparing results of numerical 
simulation with available experimental data. The exposition is organized as follows. Section 
2 deals with the kinetics of rearrangement for a temporary network of long chains. Stress- 
strain relations for a transient polymeric network are derived in Section 3. The constitutive 
equations are applied to describe uniaxial extension of a specimen in Section 4. Section 5 is 
concerned with comparison of experimental data for filled rubbery polymers with results of 
numerical simulation. Some concluding remarks are formulated in Section 6. 

2. Rearrangement of a temporary network 

The distribution of chains in a polymeric network is determined by the total number of 
active chains per unit mass, S, and the probability density of long chains composed of n 
strands, p(n). We focus on the stationary response of filled elastomers in dynamic mechanical 
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tests with small strains and assume that the quantities S and p(n) are independent of time 
(however, these parameters may, in general, depend on the amplitude of oscillations). 

Denote by X(t, r, n) the number of active chains with n strands (per unit mass) at time t 
that has last been linked to the network at instant r e [0, t]. The function X(t, r, n) entirely 
determines the current state of a polymeric network: 

X(t,t,n) = Ep(n) (3) 
equals the number of active chains (per unit mass) with n strands at time t; 
dr 



dX. 

-fr(t,T,n) 



t=T 



is the number of active chains with n strands (per unit mass) that have merged with the 
network within the interval [r, r + dr] ; 

dX f 

— (t,T,n) dr 

is the number of these chains that did not break during the interval [r,t\; 

dX, n s , 

- — (t,0,n) dt 

is the number of active chains (per unit mass) that detach from the network (for the first 
time) within the interval [t,t + dt], and 

~{t^n) dtdr 

is the number of chains (per unit mass) that has last been linked to the network within the 
interval [r, r + dr] and detach from the network (for the first time after merging) during the 
interval [t, t + dt]. 

The kinetics of reformation is determined by the relative rate of breakage for active 
chains, T(n), and by the rate of merging of dangling links with the network, 7(n). The 
relative rate of breakage is defined as the ratio of the number of active chains broken per 
unit time to the total number of active chains, 



r(n) 



1 



dX 



X(t,0,n) dt 



(t, 0, n) = - [— (t, t, n)\ g^{t, r, n) 



(4) 



The rate of breakage T may, in general, depend on the current time t, on the instant r 
when a chain has been last bridged to the network and on the guiding vector 1 (the unit 
vector whose direction coincides with the end-to-end vector R of the chain at time r). The 
neglect of these dependences means that the network is homogeneous (the rate of breakage 
for a chain is independent of the instant of its connection to the network, r, and on the 
current time, t) and isotropic (the rate of breakage is independent of the chain direction at 
the instant of merging). For a homogeneous and isotropic network, the rate of merging, 7, 
is independent of the instants t and r, as well as of the vector 1. It is defined as the number 
of dangling chains (per unit mass) bridged to the network per unit time, 



7(n) = -fr{t,T,n) 



T=t 



(5) 



Integration of differential equations (4) with the initial conditions (3) and (5) implies that 



X(t,0,n) = Sp(n)exp -T(n)t 



-^(t,T,n) = 7(n)exp[-r(n)(f - r) 



(6) 



Because the total number of active chains remains time- independent, the number of broken 
chains coincides with the number of reformed chains 

dX , .ft d 2 X . 
7(n) = -frit ,0,n) ~ I ^,r,n)dr. 

Substitution of Eq. (6) into this equality implies the balance law 
7(71) = r(n)|sp(n) exp —T(n)t +7(71) J exp —T{n)(t — t) dr\. 
The solution of this equation reads 



7(71) = ET(n)p(n). 



(7) 



Equations (6) and (7) entirely determine the current state for a temporary network of long 
chains. 
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3. Constitutive relations for a transient network 



The stress in a dangling chain is assumed to totally relax before this chain captures a new 
junction, which implies that the natural (stress-free) state of an active chain coincides with 
the deformed state of the network at the instant of their merging. The extension ratio for a 
chain bridged to the network at time r and not broken within the interval [r, t] reads^l 

X(t,r,l) = [l-C(i,r)-l]i 

where C(t,r) is the relative Cauchy deformation tensor for transition from the deformed 
state of the network at time r to the deformed state at time t and the dot stands for inner 
product. At small strains, this equality implies that 

e(t, r, 1) = 1 • £(t, t) • 1, £{t, r) = e(t) - e(r), (8) 

where e = A — 1 is the longitudinal strain for a chain, £ = |(C — I) is the relative Cauchy 
strain tensor for the network, I is the unit tensor and e(t) is the strain tensor for transition 
from the stress-free state of the network to its deformed state at time t. By analogy with the 
Kratky-Porod model, we suppose that strains in all strands coincide, which implies that the 
strain per strand equals e/n for a chain comprising n strands. Modeling strands as linear 
elastic solids with a constant rigidity /i, we find the mechanical energy of a chain, 

With reference to the conventional assumption that the excluded-volume effect and other 
multi-chain effects are screened for an individual chain by surrounding macromolecules, we 
neglect the energy of interaction between chains and calculate the mechanical energy of a 
network as the sum of the mechanical energies for individual chains. The mechanical energy 
(per unit mass) of initial chains which have not been broken within the interval [0, t] is given 
by 

^E |;*(*.0>") J/(W)dA(\), (9) 



where S is the unit sphere in the space of vectors 1 and oL4(l) is the surface element on S. 
The mechanical energy (per unit mass) of chains that merged with the network within the 
interval [r, r + dr] and are connected with the network until the current time t reads 



(t,T,n)dr J e 2 (t,r, l)dA(l). 



(10) 



An ^ 2n dr 

To calculate the integrals over S, we introduce a Cartesian coordinate frame {xi} whose 
base vectors kj are directed along the eigenvectors of the symmetrical tensor 5: 



£ = £iki + 5 2 k 2 + 5 3 k 3 , 

where 5j is the ith eigenvalue of the tensor £. The position of the unit vector 1 with respect 
to the coordinate frame is determined by the Euler angles f and 

1 = sint?(cos(/9 ki + siny? k 2 ) + cosfi k 3 . 

Substitution of these expressions into Eq. (8) results in 

e = sin 2 cos 2 (p + 5 2 sin 2 ip) + 5 3 cos 2 



Bearing in mind that dA = sin fidfidip, we obtain 



e 2 dA = / sin fid-d 



o 

47T 







2tt r , 2 

sin 2 $(5i cos 2 + 5 2 sin 2 </?) + £ 3 cos 2 $ 



/ 1 -7T r- 

= — [3(5 2 + 5 2 + 5 3 2 ) + 2(5x52 + 5i5 3 + 5 2 5 3 ) 



15 

Summing expressions (9) and (10) and using this equality, we arrive at the formula for the 
strain energy density of the network (per unit mass) 

u tt) = £ £ { X ^ °» ») f 3i ?(*(*» °)) " ^(^(^ °))" 

dU n=l L 

+ X ^ r (t,r,n)[3/ 2 (5(t,r))-4/ 2 (5(t,r))]rfr|, 

where stands for the ith principal invariant of a tensor. It follows from this equality, the 
incompressibility condition, h{£) = 0, and Eq. (8) that 
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^) = "jf E -[^,0,n)/ 2 (e d (t)) + j( _(t, r ,n)/ 2 ( ed (t) -e d (r))dr 



(11) 



where e d is the deviatoric part of the strain tensor e. The deviatoric component a d of the 
stress tensor a is determined by the conventional relationship 
dUit) 

ad{t) = P d^ry 

where p is mass density. Substitution of Eqs. (6), (7) and (11) into this equality results in 
the stress-strain relation 



a(t) = -P(t)l + ^ppE £ ^ exp(-I»t)e d (t) 



n=l 



+r(n) J exp(-r(ra)(f - r)) (e d (t) - e d (r))rfr 
where P is pressure. 



(12) 



4. Uniaxial extension of a network 

For uniaxial extension of an incompressible network along the axis x\ of the Cartesian 
coordinate frame {xi}, the strain tensor reads 

e(t) = e(f)[k!k! - i(k 2 k 2 + k 3 k 3 )], 

where e(t) is the longitudinal strain. Substitution of this equality into Eq. (12) results in 
the formula for the longitudinal stress cr(t), 
1 _ » p(n) 



a(t) = -ppE £ — e(t) - r(n) / exp(-I»(t - r))e(r)dr 

5 n=l U L 70 V 7 



(13) 



We now aim to established a correspondence between the longitudinal strain for the network, 
e(t), and the macro-strain for a specimen, 
L(t) - L(0) 



e(f) = 



L(0) 



(14) 



where L(t) is the current length of a specimen. At low concentrations of filler (the DC 
regime), we assume that the host material is homogeneous and postulate that the micro- 
strains in active chains coincide with the macro-strain in the specimen, 
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e(t) = e(t). (15) 

At high degrees of loading (the NR regime), deformation of a part of the rubbery polymer 

is assumed to be screened by the network. To describe restrictions imposed by the filler 

network on deformation of a reinforced elastomer, we (i) postulate that the network is rigid 

and (ii) characterize the level of constrains by the volume fraction of the undeformed host 

material. For uniaxial loading, this fraction is uniquely determined by the length, L Q (t), of 

a part of the sample occupied by occluded rubber. Under the hypothesis that deformation 

in the rest of the specimen is homogeneous, we find the longitudinal elongation for the 

unconstrained part of the bulk material as the ratio of the current length of this part to its 

initial length, 

_ L(t) - L (t) 
M ) ~ L(0)-L (ty 

We substitute Eq. (14) into this equality, introduce the longitudinal strain by the conven- 
tional formula, e# = A* — 1, neglect terms of the second order of smallness compared with 
e(t), and postulate that the micro-strains in long chains coincide with the macro-strain in 
the unrestricted part of the rubbery polymer. This implies that 

e(t) = ZtoF^bf W' < 16 » 

When the filler network is absent, Lo(t) = 0, and Eq. (16) turns into Eq. (15). A similar 
equality was proposed in Ref.S based on the concept of the strain-amplification factor. Un- 
like previous studies, where the coefficient on the right-hand side of Eq. (16) was expressed 
in terms of the filler concentration and the ratio of the average length to the average width 
of chain-like filler clusters by using a phenomenological equation, we do not concretize the 
shape of the function L (t) and determine it by fitting observations. As will be demon- 
strated in Section 5, this function strongly depends on the strain intensity, which implies 
that treatment of Lq as a function of the volume fraction of filler is an oversimplification. 
For a dynamic test with an amplitude A and a frequency u, we set 

e(t) = Aexp(L;t), (17) 
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where i = v / ~T- Combining Eqs. (13), (16) and (17) and introducing the new variable 
s = t — t, we find that 



E . M = *„« [i - r W /'exp(-(r W+i ,) s ) L( ^- o ^ ) 8) ^ 

where 

is the complex modulus and K (t) = |/ipSL(0)[L(0) — L (i)] _1 . 

It is natural to assume that under periodic loading, the function L (t) rapidly approaches 
its limiting value Loo(A) corresponding to the steady-state regime of oscillations. Replacing 
the function L (t) by the constant Loo (A), putting t = oo as the upper limit of integration, 
and splitting E* into the sum 

E* = E' + iL", 

where E'{uj) is the storage modulus and E"{uj) is the loss modulus, we arrive at the formulas 
^n[P(n)+^]' ^n[P(n)+^] 

with 

« - 

Given quantities N and N±, Eqs. (18) are determined by two adjustable parameters, K 
and S, and the rate of slippage from temporary junctions, T{n). To reduce the number of 
constants to be found, we set 

r(n) = T n K , (20) 

where T and /t are positive parameters. To explain the increase in the rate of breakage 
for active chains with the number of strands, we assume that chains have finite bending 
stiffness. A semiflexible chain is treated as a curvilinear rod whose micro-motion is confined 
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L(0) - L (t) 



to some tube formed by surrounding macromoleculeso. When the rod has a negligible 
bending rigidity, lateral oscillations driven by thermal fluctuations have a local character, 
which implies that slippage of a chain end from a junction is induced by random excitations 
occurring (rather sparsely) in the close vicinity of this end. For a rod with a finite bending 
rigidity, lateral oscillations induce transverse waves, whose amplitudes increase because of 
their interaction with end-points. The amplification of random transverse oscillations of an 
elastic rod in the vicinity of its end-points results in the growth of the rate of slippage of chain 
ends from temporary junction. The account for finite bending stiffness of semiflexible chains 
allows the monotonic growth of the function T(n) to be explained. It does not, however, 
result in the power-law dependence (20), which is treated as purely phenomenological. 



5. Comparison with experimental data 

We begin with fitting observations in dynamic tensile tests for a carbon black (CB) filled 
natural rubber vulcanizate with the degree of loading <p = 50 phr (parts per hundred parts 
of rubber) at various temperatures T. For a detailed description of specimens and the 
experimental procedure, see RefS. To approximate experimental data, we set N = 40, 
which is in qualitative agreement with the average number of strands used in other studies 
(for example, N = 12 was assumed for glassy styrene-acrylonitrile polymer^ and N = 20 
was employed in molecular dynamics simulations of interactions between a polymeric melt 
and filler particlesil), and N% = 500. Numerical analysis demonstrates that changes in the 
values of N and N x (in the range from iV = 30 to iV = 50 and from Ni = 100 to N x = 1000) 
weakly affect the quality of fitting. 

At any temperature T, we, first, fit the curves E'{uj) measured at the minimum amplitude 
of oscillations, A min = 0.006, and determine parameters r , k and £ (by the steepest-descent 
procedure) which ensure the best quality of matching. Afterwards, we fix the values of these 
quantities and approximate observations by using only one adjustable parameter, K, which 
is determined by the least-squares algorithm. Figures 1 demonstrate fair agreement between 
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experimental data and results of numerical simulation. 

The independence of the standard deviation of the number of strands, E, from the ampli- 
tude of oscillations, A, implies that the viscoelastic response of the filled rubber corresponds 
to the DC regime (low degrees of loading) at all temperatures T. 

The elastic modulus K is plotted versus the amplitude of oscillations A in Figure 2A 
which demonstrates that the modulus slowly decreases with the strain amplitude. Experi- 
mental data are fairly well approximated by the function 

\ogK = a -a 1 A, (21) 

where adjustable parameters a, are found by the least-squares technique. According to 
Eq. (19), below the cluster-network transition point, the modulus K is proportional to the 
number of active chains per unit mass, S. A decrease in K with A means that the growth 
of strains results in a reduction of the number of active chains because of the mechanically 
induced agitation of the breakage process. This conclusion is in agreement with observations 
for unfilled rubbery polymersil. 

The quantities £ and a Q = log if (0) (the logarithm of the elastic modulus K for the 
infinitesimally small amplitude of oscillations) are plotted in Figure 2B versus the degree 
of undercooling, AT = T — T g , where T g is the glass transition temperature. Experimental 
data are approximated by the functions 

logS = 6 -A,, logtf(O) =c + ^p (22) 

where the parameters bi and q are found by the least-squares algorithm. Equations (22) 
provide good fitting of observations. An increase in the standard deviation of the number of 
strands S with temperature is in agreement with the proposed mechanism for transformation 
of the tightly bound domains into the loosely bound regions. A decrease in the elastic 
modulus with temperature confirms the hypothesis about the elastic nature of the response of 
long chains (according to the conventional theory of entropic elasticity, the modulus linearly 
increases with temperature because of the growth in the number of available configurations 
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for long chainsc5). Equations (19) and (22) imply that the rigidity of strands, //, decreases 
with temperature. 

The rate of breakage for active chains To is plotted versus the degree of undercooling 
in Figure 2C, which demonstrates that To monotonically decreases with temperature and 
reaches its limiting value, T^, at high temperatures. Experimental data are approximated 
by the phenomenological relation 

logr = d - diAT (23) 

with adjustable parameters dj. 

At first glance, the decrease in the rate of breakage with temperature is rather surprising, 
because the growth of temperature should result in an increase in the rate of thermally acti- 
vated processes. Observations depicted in Figure 2C may be adequately explained, however, 
in the framework of the model of transient networks of semiflexible chains. The growth of 
temperature implies, on the one hand, an increase in the rate of thermal fluctuations which 
intensifies the breakage process. On the other hand, it induces a decrease in the bending 
rigidity of elastic chains, which means that transverse oscillations driven by thermal fluctua- 
tions at intermediate points of a chain are not transmitted to its ends, and, as a consequence, 
their amplitudes do not grow because of interactions of bending waves with the end-points. 

The parameter k that describes the influence of the length of an active chain on its rate 
of rearrangement is depicted in Figure 2D as a function of temperature T. The quantity k 
increases with T at relatively low temperatures and becomes practically constant at elevated 
temperatures. 

We proceed with fitting observations in dynamic tensile tests at room temperature for 
CB filled bromobutyl rubber. A detailed description of specimens and the experimental 
procedure is given in RefM To reduce the number of adjustable parameters, we fix the 
average number of strands in a chain (the same as for the CB filled natural rubber), N = 40, 
and set N\ = 80 (which implies the symmetry of the probability density p(n) with respect 
to the point N). Matching observations in a tensile test with the minimum amplitude of 
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oscillations A min = 0.0003 results in r = 0.0088 s _1 . We fix this value and approximate 
experimental data in other tests with the help of the adjustable parameters n, £ and K. 
For any degree of loading, 0, we determine the parameter k by fitting the curve E'{uj) 
measured at the minimum amplitude of strains, A min . Afterwards, we fix the value of k and 
match observations in dynamic tests with other amplitudes of oscillations by using only two 
adjustable parameters, £ and K. The quantity £ is determined using the steepest-descent 
procedure and the constant K is found by the least-squares algorithm. Figure 3 demonstrates 
fair agreement between experimental data and results of numerical simulation. 

The parameters £ and K are plotted versus the amplitude of oscillations A in Figure 4. 
Experimental data are fairly well approximated by the functions 



where the parameters Ai and Bi are determined using the least-squares algorithm. 

Two different regimes in the material behavior are revealed. At small amplitudes, A < 
A cr , £ and K weakly change with A (curves 1). The standard deviation of the number of 
strands decreases at the minimum degree of loading, = 30 phr, and remains practically 
constant for > 50 phr. The elastic modulus K increases with A at < 50 phr and decreases 
at > 50 phr. These changes may be associated with transition of bound rubber from the 
tight state to the loose state. At large amplitudes, A > A cr , £ monotonically decreases and 
K strongly increases with A (curves 2), which may be explained by disintegration of the 
filler network and release of occluded rubber. 

The parameters Aq and A\ corresponding to the post-fracture regime for the filler network 
are plotted versus the degree of loading in Figure 5A. Analogous dependences for the sub- 
fracture regime are not presented because of the insufficient number of experimental data 
for the parameter £ (the best fit of observations depicted in Figures 4C and 4D is reached 
for the homogeneous distribution of chains with various numbers of strands). Experimental 
data in Figure 5A are fairly well approximated by the functions 



log£ = A - AilogA, 



logK = B + B 1 logA, 



(24) 



A = 



Ax) + A Q1 4>, 



A 1 = A w - 



A ll( j) 



(25) 



17 



where adjustable parameters Aij are found by the least-squares method. It follows from Eqs. 
(24) and (25) that in the interval of strains under consideration, the standard deviation of the 
number of strands decreases with the filler concentration, in agreement with the proposed 
scenario (at high volume fractions of filler, aggregation of clusters into a network results in 
a decrease in the number of isolated particles, and, as a consequence, in a decrease in the 
number of chains belonging to the surface layers). 

The parameter Bo and B x for the sub-fracture regime for the filler network are depicted 
in Figure 5B versus the degree of loading 0. Experimental data are correctly predicted by 
the functions 

B = B 00 - B QX <f>, B x = B 10 - B n <f>, (26) 

where adjustable parameters By are determined by the least-squares algorithm. Because the 
quantities Bq and B\ remain practically constant in the post-fracture regime, appropriate 
figures for these parameters are omitted. 

Results presented in Figures 5A and 5B may be treated as a confirmation of the model 
because they reveal a plausible (from the physical standpoint) behavior of adjustable pa- 
rameters in the constitutive equations with changes in the filler concentration. 

The critical amplitude of oscillations A cr corresponding to fracture of the filler network 
is plotted in Figure 5C versus the degree of loading <fi. Observations are fairly well approxi- 
mated by the function 

A cr = C{<f> - (27) 

where the parameters C, (3 and cr are found by the steepest-descent procedure. Following 
Ref.i, we find that the critical loading <p cr = 41 prh corresponds to the critical volume 
fraction if) CT = 0.154, which is rather close to the theoretical percolation threshold if)* = 0.18 
for a spatial network of hard spheres@, as well as to the critical volume fraction ip* = 0.14 
for the isolator-conductor transition for CB natural rubber determined in dc conductivity 
tests@. The exponent j3 = 0.662 in Eq. (27) is also in good agreement with the theoretical 



predictions /3* = 0.70 ± 0.02, see Ref.El, and /3* = 0.72, see RefM Our estimate for the 
critical loading is confirmed by observations depicted in Figures 4A and 4B, which reveal 
that transition from the DC regime to the NR regime occurs at the degree of loading, <f), 
belonging to the interval between 30 and 50 phr. 

The exponent k is plotted versus the degree of loading in Figure 5D, which demonstrates 
that the parameter k (which characterizes the rate of breakage for long chains) remains 
constant for the NR regime. 

To assess the effect of temperature T on the distribution of long chains, we fit experimen- 
tal data for the CB filled bromobutyl rubber in dynamic tensile tests at the temperatures 
T = 25, T = 50 and T = 100 °C. The amplitude of strain oscillations A = 0.01 corre- 
sponds to the regime of rupture for the filler network. Figures 6A and 6B demonstrate fair 
agreement between observations and results of numerical analysis with N = 40, iV\ = 80 
and k = 16.0. In matching experimental data at various temperatures, the parameter Tq is 
treated as an adjustable parameter. Surprisingly, it is found that its value that ensures the 
best fit of experimental data is independent of temperature and equals T = 0.0088 s _1 . 

The standard deviation of the number of strands, E, and the elastic modulus K are 
depicted in Figures 6C and 6D versus the degree of undercooling, AT. Observations are 
approximated by the functions (22) used to match experimental data for the CB filled natural 
rubber in the DC regime. Comparison of Figure 2B with Figures 6C and 6D implies that the 
materials under consideration demonstrate different behavior for the DC and NR regimes. 
In the regime of disjoint clusters, the growth of temperature results in an increase in E and 
a decrease in K. On the contrary, in the regime of network rupture, the standard deviation 
of the number of strands decreases and the elastic modulus increases with temperature. The 
rates of changes in these parameters increase with the volume fraction of filler which may 
serve as a measure of connectivity for the network. The difference in the material response 
observed in these figures is in agreement with the proposed scenario for mechanically induced 
changes in the micro-structure of filled elastomers. 

Finally, we approximate observations for styrene-butadiene rubber filled by carbon black 
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and polymeric particles in dynamic shear tests with the strain amplitude A = 0.002 at 
the temperature T = —30 °C. A detailed description of specimens and the experimental 
procedure is found in Ref.i. Our aim is to assess the influence of the filler material on 
adjustable parameters in the constitutive equations. 

Figures 7A and 7B demonstrate fair agreement between experimental data and results of 
numerical simulation with N = 40, Ni = 80, Tq = 0.0012 s _1 and k = 4.33. The quantities 
£ and K are depicted versus the degree of loading in Figures 7C and 7D. Observations 
are fairly well approximated by the functions 

log S = C + d0, log K = D - D 1( j>, (28) 

where the constants C, and Di are determined by the least-squares technique. Figure 7 
reveals transition from the DC regime of micro-structural changes at small concentrations 
of filler to the NR regime at large degrees of loading. The critical concentration reads 
4> CT ~ 40 phr for CB particles and <j) cr ~ 30 phr for polymeric particles. For the DC regime, 
the standard deviation of the number of strands, S, slowly changes with <p (increases for 
the polymeric filler and decreases for the CB filler, in agreement with data depicted in 
Figure 4A), whereas the elastic modulus K increases with for both kinds of filler (in 
agreement with observations presented in Figure 4A). For the NR regime, the quantity S 
increases and the modulus K decreases with <fi for both kinds of filler. It is worth noting 
that the rates of changes in £ and K with the degree of loading are similar for CB and 
polymeric particles. This means that in the regime of disintegration of the filler network, the 
viscoelastic response of filled elastomers is mainly determined by the topological properties 
of the network and it rather weakly depends on physical properties of the reinforcement. 

6. Concluding remarks 

Constitutive equations have been derived for the mechanical response of particle-reinforced 
elastomers at small strains. The model is applied to describe the Payne effect: changes 
in the dynamic moduli of filled rubbery polymers with an increase in the amplitude of 
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strain-oscillations. The viscoelastic behavior of reinforced rubbers demonstrates two differ- 
ent regimes associated with (i) transformations of bound rubber in the close vicinity of filler 
particles and their clusters and (ii) fracture of the filler network and its disintegration into 
disjoint clusters of particles. The critical concentration of filler corresponding to the transi- 
tion from the DC regime to the NR regime is in good agreement with theoretical predictions 
by the percolation theory, as well as with experimental data for dc conductivity. 
The following conclusions are drawn: 

1. Constitutive equations correctly predict observations in dynamic mechanical tests on 
filled elastomers at various temperatures, filler concentrations, amplitudes and fre- 
quencies of oscillations. 

2. Changes in conditions of mechanical experiments result in changes in adjustable pa- 
rameters of the model which are adequately described within the concepts of bound 
and occluded rubber and which are in agreement with the proposed scenario for tran- 
sition from the DC regime to the NR regime of deformation. 

The work is confined to the effects of temperature, filler concentration and amplitude and 
frequency of strain oscillations on the viscoelastic response of filled rubbers. Some important 
questions, however, remain beyond the scope of this study. They include, in particular, (i) 
the influence of the shape of particles, the distribution of their sizesil, and the presence of 
coupling agents0 on the rheology of reinforced elastomers, (ii) changes in the glass transition 
temperature of filled rubbery polymers caused by polymer-filler interactions, as well as (iii) 
the effect of mechanical loading on electrical conductivity of reinforced elastomersl3i^l. These 
issues will be the subject of a subsequent publication. 
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FIGURES 

Fig. 1. The storage modulus E' MPa versus the frequency of oscillations oj Hz for CB filled 
natural rubber in a tensile dynamic test with the strain amplitude A at a temperature T K. 
Circles: experimental dataS. Solid lines: results of numerical simulation. A: T = 253; B: T = 296; 
C: T = 333; D: T = 373. Curve 1: A = 0.006; curve 2: A = 0.011; curve 3: A = 0.028; curve 4: 
A = 0.056 

Fig. 2. A: The elastic modulus K versus the amplitude of oscillations A at a temperature T K. 
Solid lines: approximation of the experimental data by Eq. (21). Curve 1: T = 253, ao = 7.2672, 
oi = 5.3129; curve 2: T = 296, a = 4.6737, oi = 3.4019; curve 3: T = 333, a = 3.0929, 
a t = 2.5039; curve 4: T = 373, a = 2.8935, ai = 1.9944. B: The standard deviation of the 
number of strands S (unfilled circles) and the initial modulus K(0) (filled circles) versus the degree 
of undercooling AT K. Solid lines: approximation of the experimental data by Eq. (22). Curve 1: 
bo = 2.0465, bi = 16.246; curve 2: Co = 1.7553, c\ = 185.33. C: The rate of breakage Tq s _1 versus 
the degree of undercooling AT K. Solid line: approximation of the experimental data by Eq. (23) 
with d = -0.9827 and di = 9.3009 x 10~ 2 . D: The coefficient K versus the degree of undercooling 
AT K. Solid line: the average value k = 5.7. Symbols: treatment of observations for CB filled 
natural rubbe 

Fig. 3. The storage modulus E' MPa versus the frequency of oscillations uo Hz for CB filled 
bromobutyl rubber in a tensile dynamic test with the strain amplitude A at room temperature. 
A: 4> = 30; B: = 50; C: <f) = 70; D: 4> = 100. Circles: experimental data!. Solid lines: results of 
numerical simulation. Curve 1: A = 0.0003; curve 2: A = 0.0011; curve 3: A = 0.0033; curve 4: 
A = 0.0048; curve 5: A = 0.0077; curve 6: A = 0.0122 
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Fig. 4. The standard deviation of the number of strands S and the elastic modulus K MPa 
versus the amplitude of oscillations A for CB filled bromobutyl rubber. Circles: treatment of 
observationsi. Solid lines: approximation of experimental data by Eq. (24). A: <p = 30, Curves 1: 
A) = 0.6474, Ay = 0.1208; B = 10.0943, By = 1.5797; B: = 50, Curves 1: A = 0.9702, 
Ay = 0.0084, B = 8.4060, By = 0.7724, curves 2: A = -0.1828, Ay = 0.4658, B = 42.2894, 
By = 14.5677; C: <p = 70, Curve 1: B = 3.7590, By = -0.0458, curves 2: A = -0.0706, 
Ay = 0.3871, B = 44.4787, By = 14.4917; D: <f> = 100; Curve 1: B = 2.8452, By = -0.7079; 
curves 2: A = 0.1749, Ay = 0.2567; B = 39.7411, By = 11.5563 

Fig. 5. A: The parameters Ao (unfilled circles) and Ay (filled circles) versus the degree of loading 
(p phr. Solid lines: approximation of the experimental data by Eq. (25). Curve 1: Aqo = —0.5568, 
A y = 7.2353 x 10~ 3 ; curve 2: A 1Q = 0.6775, Ayy = 4.1950 x 10" 3 . B: The parameters B (unfilled 
circles) and By (filled circles) versus (p phr. Solid lines: approximation of the experimental data 
by Eq. (26). Curve 1: B 00 = 13.3942, B 0l = 0.1059; curve 2: B 10 = 2.4530, Byy = 0.0329. C: The 
critical amplitude of oscillations Ae cr versus <j) phr. Solid line: approximation of the experimental 
data by Eq. (27) with C = 1.4711 x 10" 2 , (5 = 0.6616 and 4> cv = 40.99 phr. D: The parameter k 
versus <j) phr. Solid line: the limiting value k = 16.0. Symbols: treatment of observations for CB 
filled bromobutyl rubber at room temperature! 

Fig. 6. A and B: The storage modulus E' MPa versus the frequency of oscillations u Hz for 
CB filled bromobutyl rubber in a tensile dynamic test with the strain amplitude A = 0.01 at a 
temperature T °C. Symbols: experimental data!. Solid lines: results of numerical simulation. A: 
(j) = 60 phr, B: (p = 100 phr. Curves 1: T = 25; curves 2: T = 50; curves 3: T = 100. C: The 
standard deviation of the number of strands £ versus the degree of undercooling AT K. Symbols: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. (22). 
Curve 1: <p = 60 phr, b = 0.6956, by = -8.8493; curve 2: <p = 100 phr, b = 0.5400, by = -18.337. 
D: The elastic modulus K MPa versus the degree of undercooling AT K. Symbols: treatment of 
observations. Solid lines: approximation of the experimental data by Eq. (22). Curve 1: <p = 60 
phr, c = 14.084, cy = -236.96; curve 2: cp = 100 phr, c = 24.086, cy = -911.59 
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Fig. 7. A and B: The storage modulus E' MPa versus the frequency of oscillations u rad/s for 
filled styrene-butadiene rubber (SBR-1500) with a degree of loading (ft phr in a tensile dynamic test 
with the strain amplitude A = 0.002 at the temperature T = —30 °C. A: CB filled elastomer; B: 
polymeric filled elastomer. Circles: experimental datai. Solid lines: results of numerical simulation. 
Curve 1: (ft = 0; curve 2: <ft = 25; curve 3: (ft = 50; curve 4: (ft = 75; curve 5: (ft = 100. C: The 
standard deviation of the number of strands £ versus the degree of loading (ft phr. Solid lines: 
approximation of the experimental data by Eq. (28). Curve la: Cq = 1.0864, C\ = —0.0013; 
curve 2a: C = 0.7832, C x = 0.0063; curve lb: C = 1.0864, C x = -0.0010; curve 2b: C = 0.9003, 
C\ = 0.0071. D: The elastic modulus K MPa versus (ft phr. Solid lines: approximation of the 
experimental data by Eq. (28). Curve la: D = 3.8947, D x = -0.0211; curve 2a: D Q = 4.8175, 
D x = 0.0064; curve lb: D = 3.8947, D x = -0.0097; curve 2b: D = 4.2435, D 1 = 0.0055. 
Symbols: treatment of observations!!. Unfilled circles: CB filler; filled circles: polymeric filler 
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